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Abstract 

Dynamics of sine-Gordon kinks in the presence of rapidly varying periodic 
perturbations of different physical origins is described analytically and nu- 
merically. The analytical approach is based on asymptotic expansions, and 
it allows to derive, in a rigorous way, an effective nonlinear equation for the 
slowly varying field component in any order of the asymptotic procedure as 
expansions in the small parameter u; -1 , uj being the frequency of the rapidly 
varying ac driving force. Three physically important examples of such a dy- 
namics, i.e., kinks driven by a direct or parametric ac force, and kinks on 
rotating and oscillating background, are analysed in detail. It is shown that 
in the main order of the asymptotic procedure the effective equation for the 
slowly varying field component is a renormalized sine-Gordon equation in the 
case of the direct driving force or rotating (but phase-locked to an external ac 
force) background, and it is the double sine-Gordon equation for the paramet- 
ric driving force. The properties of the kinks described by the renormalized 
nonlinear equations are analysed, and it is demonstrated analytically and 
numerically which kinds of physical phenomena may be expected in dealing 
with the renormalized, rather than the unrenormalized, nonlinear dynamics. 
In particular, we predict several qualitatively new effects which include, e.g., 
the perturbation-induced internal oscillations of the 2n— kink in a parametri- 
cally driven sine-Gordon model, and generation of kink motion by a pure ac 
driving force on a rotating background. 
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I. INTRODUCTION 



As is well known, the effect of rapidly varying perturbations on the dynamics of non- 
linear systems may lead to a drastic change of the system behavior in the sense of the 
averaged dynamics. In particular, large-amplitude parametric perturbations may give rise 
to a stabilization of certain types of dynamical regimes. A typical and famous example is 
a stabilization of a reverse pendulum by parametric forced oscillations of its pivot Q (see 
also the recent paper and references therein), and a similar effect may be also achieved 
by applying a direct ac driving force of large amplitude ||. Such a dynamical stabilization 
has its analog in nonlinear systems with distributed parameters supporting, in particular, 
novel types of kink solitons However, the method which is usually used to derive an 

averaged equation describing the system dynamics in the presence of rapidly varying pertur- 
bations is not rigorous and, as a matter of fact, it is not well justified. Such a method, even 
being very clear from the physical point of view, uses a splitting of slow and fast variables 
and subsequent averaging which is based, in fact, on solutions of a linearized equation for 
fast variations where the slowly varying coefficients are assumed to be constant [l],f|,f| . The 
procedure of such a linearization assumes that the amplitude of the forced (rapidly varying) 
oscillations is small, and this is certainly valid for parametrically forced oscillations far from 
the parametric resonance. For direct ac perturbations, the forced oscillations may become 
large. To describe the dynamics in an approximate way, the so-called "rotating-wave ap- 
proximation" was used without detailed mathematical justification It is necessary to 
note that the derivation of an effective averaged equation for the slowly varying field com- 
ponent is an important (and nontrivial) step of the analysis of such systems, and in many 
of the cases the corresponding equation determines the leading physical effects observed in 
the presence of rapidly varying perturbations. In all the cases it is necessary to justify the 
averaging procedure as well as to estimate the influence of the higher-order contributions. 
Unfortunately, the latter are beyond the usual averaging methods. However, as we show in 
the present paper, a rigorous analysis of the effect of rapidly varying periodic perturbations 
on nonlinear dynamics of ac driven damped systems can be performed in a straightforward 
way to describe the averaged dynamics with any accuracy. 

The purpose of this paper is to present the basic steps of the method mentioned above 
and, selecting the sine-Gordon (SG) model as a particular example, to describe the dynamics 
of kinks in the presence of rapidly varying driving forces of very different physical origins. 
Considering the external ac driving force to be rapidly oscillating, we apply an asymptotic 
procedure based on a Fourier series where the coefficients are assumed to be slowly varying 
functions on the time scale a; -1 , u being the frequency of the rapidly varying ac force which 
is assumed to be large. The basic idea to split fast and slow variables is not new, and the 
well-known example is, as mentioned, a stabilization of the reverse pendulum by oscillations 
of its suspension point. However, our analytical method to derive an effective equation for 
the slowly varying field component is novel, and this method allows us to calculate, in a 
self-consistent way, all the corrections using solely asymptotic expansions rather than direct 
averaging in fast oscillations. It is clear that the applicability of the method itself is much 
wider than the particular examples covered by the present paper. 

The paper is organized as follows. In Sec. II we consider the case of a direct ac driving 
force demonstrating the basic steps of our analytical approach in detail. The main result 
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of such an analysis is the so-called "averaged" equation, i.e., that describing slowly varying 
system dynamics. In the case of the direct driving force this equation is shown to be a 
renormalized SG equation. Section III presents the case of a parametric driving force where 
the final equation describing the slowly varying dynamics is the double SG equation which, 
as we show, may display new features in the averaged kink dynamics, e.g., oscillations of 
the excited internal mode of the kink, which is absent in the standard SG model. The case 
of kink stabilization on a rotating background by applying a rapidly oscillating ac force is 
discussed in Sec. IV. There we analyse an effect of an induced dc force on the kink motion, 
as well as showing numerically that the main conclusions of our analysis may be easily 
extended to cover multi-soliton dynamics. Finally, Sec. V concludes the paper. 



II. DIRECT DRIVING FORCE 

A. Asymptotic expansions 

As the first example, we consider the case of the direct driving force in the SG model 
when the system dynamics is described by the driven damped SG equation for the field 
variable <f)(x, t), 

^-^ + W = / - 7 - + e cosM), (1) 

where / is a constant contribution of the driving force, 7 is the damping coefficient, and 
the amplitude e of the driving force may be large (in fact, up to the values of order of 
oj 2 ). The standard physical application of the model flU) is to describe the fluxon dynamics 
in long Josephson junctions (see, e.g., |§), so that / and ecos(u;t) are the constant and 
varying components of the bias current applied to the junction. In the subsequent analysis 
we consider the direct driving force (~ e) as rapidly oscillating, i.e., the frequency u is 
assumed to be large in comparison with the frequency gap (= 1) of the linear spectrum 
band. Our purpose is to derive an "averaged" nonlinear equation to describe the slowly 
varying dynamics of the SG field. 

In order to derive an averaged equation of motion, we note first that in the case of very 
different time scales the SG field (f) may be decomposed into a sum of slowly and rapidly 
varying parts, i.e., 

= $ + C (2) 

The function ( stands for fast oscillations around the slowly varying envelope function <3>, 
and the mean value of £ during an oscillation period is assumed to be zero so that < >= $. 
Our goal is to derive an effective equation for the function $. The standard way to do that 
is to substitute Eq. (||]) into Eq. (|l|) and to split Eq. ([]]) into two equations for slow and fast 
variables, making an averaging to obtain the equation for the slowly varying field component 
(see, e.g., [§,f|). However, such an approach must be properly justified for the case when the 
fast oscillations are not small as it is for the direct driving force considered here, and in a 
similar problem it was proposed |1 to use the so-called "rotating wave approximation" to find 
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the rapidly oscillating field component. All these approaches, although quite satisfactory 
for the first-order approximation (see, e.g., 0||[7]]), do not allow to make the next-order 
expansions to calculate higher-order corrections, and thus they cannot be rigorously justified. 
Nevertheless, as we show in the present paper, a rigorous approach may indeed be proposed 
to obtain an effective equation for the slowly varying field component $ with any accuracy. 

The basis of our asymptotic procedure is a Fourier series expansion with slowly varying 
coefficients. We look for rapidly oscillating component ( in the form 

( = A cos(ut) + B sin(cjt) + C cos(2cji) + D sin(2cjt) + . . . , (3) 

where the coefficients A, B, . . . are assumed to be slowly varying on the time scale ~ uj^ 1 . 
Substituting the expressions (|j), @ into Eq. (P, we note that the effective coupling between 
different harmonics of the expansion (|2|), (|3|) is produced by the nonlinear term sin0, which 
generates the following Fourier expansion, 

sin (ft = sin $ [a + a>i cos{ut) + a 2 sin(u;t) + a 3 cos(2a;t) + a 4 sin(2co>;£) + . . .] 

+ cos $ [/3q + /?i cos(cjt) + /3 2 sm(ut) + (3 3 cos(2cut) + /3 4 sin(2c^) + ...], (4) 

where 

a = J (A) (l - + \b 2 J 2 (A) + ..., a 1 = -CJ X {A) + ..., 

a 2 = -DJ 1 (A) + ..., a 3 = ]B 2 J (A) + ..., a 4 = -BJ^A) + . . . , (5) 



P = -CJ 2 (A) + ..., fa = 2J X {A) + . . . , p 2 = BJ (A) + BCJ 1 (A) + ... 

(3 3 = CJo(A) + . . . , & = DJ (A) + ... 



(6) 



and Jq, J\, etc., are Bessel functions. Collecting now the coefficients in front of the different 
harmonics, we obtain an infinite set of coupled nonlinear equations, 



d 2 $ d 2 $ 



dt 2 dx 2 



+ sin$ 



1 



4 



Jo(A) [1--B 1 )+ B 2 J 2 (A) + ... 
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+ cos *[-CJ 2 {A) + ...]=/ -7^, 



(7) 



-ufA + 2u— + 



dB d 2 A\ d 2 A 



dt dt 2 



dx 2 



+ cos$[2Ji(A) + . 



'OA 



+ sin $[-C Ji (A) + . . .] + 7 [ — + coB 



(8) 



-uj z B - 2cu— + 



dA d 2 B\ d 2 B 



dt dt 2 



dx 2 



+ cos $[B J (A) + BC Ji(A) + . . .] 
+ sin $[-£) Ji(A) + . . .] + 7 ( — - coAj = 0, 



(9) 
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-4urC + Auj— + 



3D 8 2 C\ d 2 C 



dt at 2 



dx 2 



+ cos $[CJ (A) + 



+ sin$ 



\b 2 ua) + ... 



(10) 



dC d 2 D s 



-Auj z D - Au— + 



dt dt 2 



d 2 D 

dx 2 



+ cos $[DJ (A) + 



+ sin$[-BJi(A) + ...]+ T f — 



0. 



(11) 



and similar equations for the coefficients in front of the higher-order harmonics. To proceed 
further, we note that different terms in Eqs. (§) to (|TT| ) are not equivalent provided u is a 
large parameter. Indeed, if we assume the amplitude e large as well (otherwise, the dynamics 
of the system is rather trivial because the effect of small-amplitude but rapidly oscillating 
force is negligible), let us say up to the order of uj 2 , the large term — uj 2 A in Eq. (§) may 
be compensated only by the term e from the right-hand side of Eq. (^). Thus, assuming 
e ~ uj 2 we find the first term of the asymptotic expansion A rs —e/uo 2 . On the other hand, 
the right-hand side of Eq. @ is zero, and the large term —uj 2 B may be compensated only 
by a contribution from the other terms ~ A, thus giving the first term of the expansion for 
the coefficient B, viz., B « 'ye/uj 3 . Such a simple reasoning may be effectively applied to 
other coefficients as well as to other corrections of the asymptotic expansion. As a matter 
of fact, to generalize and simplify the procedure of calculation of the expansion coefficients, 
we look for the coefficients A, B, . . . in the form of the power series in the small parameter 
uj^ 1 as follows 



A , 02 D & 1 , & 2 n Ci C 2 

A = ai + — + ..., 5 = h — + . . . , C = — + — r + D 



UJ UJ° UJ"* UJ" UJ° UJ 



4 • 4 • ••• • (12) 



Substituting Eq. ([TJ) into Eqs. (@) to (p]) and equating the terms of the same orders in 
the small parameter uj~ 1 , we find 



ai 



-5, 



(13) 



a 2 = 7&i + 2 cos $ Ji(ai), 



(14) 



6. = -2 



<9a 2 



6i = -701 



7a 2 + 6i cos$[J (ai) + 2J 2 (ai)], 



(15) 
(16) 



ci 



4^ + -6?J (ai)sin$ + -6?J 2 (ai) + 2 7 d! 



(17) 
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d\ = — ~b\ J\{a\) sin 



(18) 



and so on. In Eq. ([y|) the parameter 5 = e/u> 2 is assumed to be of order of 0(1), but all the 
results are valid also for the case 8 <C 1. The expansions (|T2| ) allow to find the coefficients 
A,B, . . . in each order of co> _1 , and all the corrections are determined by algebraic relations 
rather than additional differential equations. For example, 02 is determined by Eq. ( p"4|) 
through h\ which, in its turn, may be found from Eq. (|T5| ) as a function of a\, i.e., through the 
slowly varying part $, and so on. This statement is valid for all coefficients of the asymptotic 
expansion: The coefficients are found through algebraic relations involving lower-order terms 
of the asymptotic expansions and their derivatives, and one does not need to find solutions 
of additional differential equations. 



B. Renormalized equation 



Applying the expansions (12) to Eq. (|7|), we may find the equation for the slowly varying 
field component $ with any accuracy in the small parameter ui , e.g., 



<9 2 $ d 2 $ 



dt 2 dx 2 



+ sin$ 



Jo(ai) + -=■ f-7& 2 J (ai) - o 2 Ji(ai) + ^6 2 J 2 (ai) ) + . . . 



4 



4 



+ cos$ 



Cl 



— J 2 (ai) + ... 



uJ 1 



<9$ 

f-y w (19) 



Thus, from the asymptotic procedure described above it is quite obvious how to calculate 
the corrections of the first, second, and subsequent orders and to find the averaged equation 
with any required accuracy. 

In the first-order approximation in u~ x only the term J (ai) sin $ contributes, so that 
Eq. fllD yields 



<9 2 $ <9 2 $ 



dt 2 



dx 2 



0J< 



sin $ = / — 7 



<9$ 
~dt' 



(20) 



Equation (|20| ) takes into account an effective contribution of the rapidly varying force to 
the average nonlinear dynamics and this contribution might become large for 5 = 0(1), i.e. 
when e ~ uo 2 . Thus, the dynamics of the SG model with a rapidly varying direct driving 
force may be described by a renormalized SG equation (|20"D up to the terms of order of t/uo 2 . 

The results obtained above may immediately be applied to describe the renormalized 
dynamics of kinks in the presence of the rapidly varying ac force. In fact, Eq. (p0|) is the 
dc driven damped SG equation with a renormalized coefficient in front of the term ~ sin$. 
This simply means that we can apply all the results known for the standard SG equation 
(see, e.g., P,|i~0f) making only a renormalization of the kink's width. For example, the kink 
solution of Eq. (^) at 7 = / = has the form 



$(x, t) = 4cr tan 1 exp 



x-Vt 



(21) 
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where a = ±1 is the kink's polarity and Iq = [Jo^e/oj 2 )]^ 1 / 2 is the kink's width at rest. The 
motion of the kink in the presence of small dc force / and damping (~ 7) is characterized 
by the steady-state velocity 

V. = -', ^(ilW^V (22) 



•A + sr 

In the theory of long Josephson junctions the kink's velocity is connected with the voltage 
across the junction, < <p t >, where < . . . > stands for the averaging in time, so that the 
result ( |22"D for the steady-state kink velocity gives the so-called zero-field steps in the current- 
voltage (IV) characteristics of a long junction. As follows from Eq. (|22|) , the renormalization 
of the parameter g leads to a change of the kink's velocity V*(f) and this, therefore, changes 
the slopes of the voltage steps by the effect of the ac driving force. 



III. PARAMETRIC DRIVING FORCE 

Let us consider now a parametric driving force applied to the SG system, with the 
main purpose to demonstrate that such a case is very different from that analysed above. 
The qualitative difference between the effects produced by direct and parametric (rapidly 
oscillating) forces is the following: A sufficient change of the system "averaged" dynamics 
due to a rapidly oscillating direct force may be observed for amplitudes e ~ uj 2 whereas in 
the case of a parametric force, similar effects may be already observed for smaller amplitude, 
i.e., in fact for e ~ uj. To prove this statement and to show how our asymptotic method 
works for the case of the parametric force, we consider the parametrically perturbed SG 
equation in the form 

d 2 <j) d 2 (f) d<j) . 

where / and 7 have the same sense as above, but this time e is the amplitude of the 
parametric force. Various applications of the model ( 231 ) were discussed in the review paper 
(see also Refs. Q). 



We assume that the parametric force is rapidly oscillating, i.e., the frequency uj is large. 
As above, we look for a solution of Eq. (^3j) in the form of asymptotic expansion 

= $ + A cos(ujt) + B sin(cut) + C cos(2ujt) + D sin(2ujt) + . . . , (24) 

where the functions $, A, B, . . . are assumed to be slowly varying on the time scale ~ uj^ 1 . 
The function <3> in Eq. ( p4|) determines, in fact, the evolution of the averaged field component 
because <</>>=$, where the brackets < . . . > stand for the averaging in fast oscillations. 
Substituting the expression flU) into Eq. fl23|) and collecting, as in the case of the direct 



driving force, all the coefficients in front of the different harmonics, we again obtain an 
infinite set of coupled nonlinear equations. The subsequent (and very important) step of 
such an analysis is to find the form of the asymptotic expansions for the coefficients A, B, . . .. 
In the present case it is easy to check that the expansions (|T2"D do not give a closed asymptotic 
procedure, and in the case e ~ uj 2 the driving force from Eq. (E3|) contributes to all the 
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harmonics, so that contributions of the other harmonics become large as well. Comparing, 
as in the previous case, different terms of the equations for the coefficients A, B, . . ., we may 
easily check that the asymptotic procedure may be effectively formulated for smaller (but 
not small) amplitudes, i.e., when e ~ u, and, as above, it gives all the corrections to the 
averaged nonlinear dynamics in a rigorous way. Thus, we take the asymptotic expansions 
in the form 

B = VV„„ C=V..., + (25) 

OJ 1 CU 4 U) 6 U) b UJ^ L) b 



Using the power series fl25[) , we obtain the "averaged" equation in the form 
8 2 <t> <9 2 $ / la 2 \ . ^ , <9$ e /ai a 2 



dt 2 




sin$ = / _ 7 ^ + £eo S$ (- + 3 + ...), (26) 



The expansions (|2q) allow to find the coefficients of the asymptotic expansions in each order 
in the small parameter uj~ 1 , and all the corrections are determined, as above, by algebraic 
relations. 

In the first-order approximation only the term ~ ea\ contributes to Eq. fl26|) . From the 
asymptotic expansions it follows that 



ai = -esin$, (27) 



and Eq. fl26|) yields 



l + -A 2 cos$ sin$ = /-7— , (28) 



dt 2 dx 2 V 2 / J ' dt' 

where A = e/u. Equation (|28|) takes into account an effective contribution of the rapidly 
varying parametric force to the slowly varying nonlinear dynamics in the lowest order, and 
all the corrections coming from the approximation of the next order are proportional to the 
small parameter ~ u~~ 2 . However, even the lowest-order contribution might become large 
for A = 0(1), i.e., when e ~ to. 

Thus, the "averaged" dynamics of the SG model with a rapidly varying parametric force 
is described by the double SG equation (p8|). As a matter of fact, the double SG equation is 
rather well studied (see, e.g., Refs. [fl2| , [T~3|1 and references therein) and properties of its kink 
solutions are known as well. In particular, the kink solution of Eq. (53) at / = 7 = may 



be written in the form [12 



<fr(x,t) = 2tan" 



cosech 



1 + A 2 /2 




(29) 



and this solution may be treated as two coupled 7r-kinks. In Fig. 1 we show the results of 
numerical simulations of the parametrically driven SG system, Eq. (23). In all the cases 
analysed in the present paper we have integrated the driven damped SG equation on a 
spatial interval of length L, with periodic boundary conditions. As is seen from Fig. 1, 
the sech-type shape of the 2n kink corresponding to the standard (unperturbed) SG system 
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is modified, and the function <f> x displays a two-peaked profile which, as a matter of fact, 
is one of the main features of the kink solution (29). Increasing A 2 one may observe, in 
accordance with Eq. (p9[), that the function has an evident shape of two 7r-kinks separated 
by a distance ~ A 2 . 

As has been shown in Ref. [|J, it— kinks themselves may exist in the parametrically 
driven SG chain provided the condition A 2 > 2 is satisfied. This condition simply means 
that the effective "averaged" potential for the slowly varying field component $ exhibits a 
local minimum at $ = 7r so that this stationary state becomes stable. 

The appearence of new features in the slowly varying ("averaged") system dynamics of 
the SG system for A 2 > 2 is similar to the phenomenon of the parametric stabilization of the 
reverse pendulum in the well known Kapitza problem |],|2|]. However, in the problem under 
consideration some novel features in the nonlinear dynamics of the parametrically driven 
SG system may be really observed for any value of the effective parameter A 2 . Indeed, 
as is known from the theory of the double SG equation |T5[ , at any value of A 2 the kink 
(f29p possesses the so-called internal ( "shape" ) mode which describes variations of the kink's 
width. This internal mode is absent for the standard SG kink, and the mode frequency Q 2 sol 
splits at any A 2 ^ from the gap frequency of the linear spectrum. For A 2 > 2 the kink's 
internal mode may be described as relative oscillations of the it— kinks of which the 2tt— kink 
consists. However, this mode does, in fact, exist at any value of the effective parameter A 2 , 
and it may be observed as periodic variations of the kink's shape. 

We have measured numerically the frequency of the shape oscillations of the 27r-kink (|29"D 
directly solving Eq. (|23|) and also using the averaged equation (|8|). The numerical results 
are shown in Fig. 2 for selected values of the external frequency, u = 50 and u = 100. 
For relatively small A 2 (i.e., e in Fig. 2), when higher-order corrections to Eq. fl28|) are 
negligible, a perfect agreement between the results for the parametrically driven SG model 
(|3|) and those for the averaged equation (|^) are clearly observed, justifying the validity of 
our asymptotic procedure. 



IV. KINKS ON ROTATING AND OSCILLATING BACKGROUNDS 

As was mentioned in Ref. [0], the other physically important case when a rapidly varying 
ac force may change drastically the kink dynamics is the case of a rotating and oscillating 
background. We should note, however, that if one considers relatively small system lengths 
L, even a relatively weak driving force may lead to complicated dynamics involving coexisting 



states of bunched kinks and nontrivial background states These latter effects are 

probably caused by the influence of nonzero boundary conditions which may "help to lock" 
kink-like states on rotating backgrounds. Here we are interested in the dynamics of the 
long SG systems (i.e. the kink's length is much smaller than the system length L) when the 
high-frequency force phase-locks the SG field in an oscillating and rotating state and thereby 
creates a mechanizm (an effective gravitation field) for supporting kink solitons. However, 
we should note that the theory presented below can not be applied to formally infinite SG 
system because for the case of the continum linear spectrum the applied ac force may create 
resonances making the system dynamics much more complicated and even chaotic. In fact, 
we need finite-width (but of large L) systems in order to avoid linear resonances if the 
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frequency of the external ac force is selected in a gap between the nearest eigenfrequencies. 

To describe the effect of the kink phase-locking on a rotating background analytically in 
a rigorous way, we consider the perturbed SG equation ([D) assuming / > 1, in which case 
the ground state of the SG chain is not stable and the chain rotates with the frequency fi 
so that <p m fit. Applying the high-frequency ac force ~ e we are interested in the slowly 
varying phase-locked system dynamics on such a rotating background. Accordingly, we look 
for a solution of Eq. (|I|) in the form 

(J) = $ + fit + f , (30) 

where £ is the rapidly varying part oscillating with the large frequency u of the external 
ac force, $ is the slowly varying (long time scale) part, and fi is the average frequency of 
rotation for the background field, which we assume to be phase-locked to the external ac 
field, i.e. fi = ±ku, k being integer. Looking for the rapidly oscillating part £ in the form 
of a Fourier series with slowly varying coefficients, 

£ = A cos(ut) + B sra(wt) + . . . , (31) 

we obtain the following equations for the averaged field component $ and the expansion 
coefficients A,B, . . ., 

w-w +MA),6n * =, - ,ku - , w (32) 

/ f) R r) 2 A\ rfi A ( r) A \ 

- 2u— + — - _ + [J k+1 (A) - J k _ x {A)\ cos $ + 7 — + c5 = e , (33) 



~" 2B + 2 ^ + ^)-^- ^» iA > + J ^ 8il1 * + ^ ( St ~ uA ) = °' (34) 

and so on. Unlike the case considered in Sec. II, in the present problem there are two 
rapidly oscillating contributions with the frequencies u and ku, so that the final equations 
(j32|) to ( j34|) for the slowly varying coefficients differ from the corresponding equations (0) 
to (H of Sec. II. To take into account the second term in the right-hand side of Eq. (|32| ) in 
a self-consistent way, we also assume the dissipation to be rather small, jlu ~ 1, which is a 
typical case for Josephson junctions. 

Making now asymptotic expansions similar to the case of the direct ac force considered 
above, i.e., 

A = a 1 + ^ + ..., S = ^ + ..., (35) 



we find the following relations [cf. Eqs. (0), (|15|)] 



ai = — e — = -5, bi = -7<2i, (36) 



10 



which allow to obtain the effective equation for the slowly varying system dynamics by just 
combining Eqs. (|36l), (|35|) and fl32[). The final equation is just Eq. fl32|) with A = —e/u> 2 , 
which describes the kink dynamics on the background rotating with the frequency Q = ±fca>. 
It is important to note that the resulting (effective) dc force in the "averaged" nonlinear 
equation (|32|) is represented by the term f—jkoj but not / itself, i. e., the kink on the rotating 
and oscillating background may move even in the absence of the constant contribution to 
the bias current, / = 0. Figure 3 shows the results of the numerical calculation of the first 
Shapiro step (k = 1) of a long Josephson junction described by the model (^) with the 
parameters 7 = 0.1, e = 12.5, and Q = uj = 2.5. As is clearly seen from that figure, the 
steps cross the zero current axis, displaying the property mentioned above; constant-voltage 
zero-crossing steps are of considerable practical interest for voltage-standard applications 
of Josephson junctions (see, e.g., JT3]| and references therein). If we select / = 0.25, the 



effective force acting on the kink vanishes, and the kink is observed at rest (see Fig. 4). 
This result is in excellent agreement with the averaged SG equation (|3"2]), where the effective 
force is found to be / — ^kuj. Increasing the value of the bias current /, we create an effective 
force acting on the kink and it moves to the left (see Fig. 5). It is very interesting to note 
that the kink motion is even possible in the absence of the constant component of the bias 
current, i.e., at / = 0. Figure 6 shows this case, when the effective force — , ykuj generates 
the motion of the kink in the direction opposite to that shown in Fig. 5. 

As a matter of fact, our approach and the resulting equation for the slowly varying field 
component $ do not specify exactly the type of nonlinear solutions we deal with. This 
means that the consideration described above may be effectively applied to other problems, 
much more general than a single-kink propagation. One of the important generalizations of 
this approach is to treat multi-soliton (multi-kink) problems. In particular, for the nonlinear 
dynamics on oscillating and rotating background the effects similar to those described for a 
single kink may be also observed for the case of two, three, and more kinks. In particular, 
Figures 7 and 8 present results of our numerical simulations of the same effects as for the 
single kink, but for the case of two kinks in the directly driven SG model ([J). As may 
be noted from these figures, the large-amplitude direct ac force generates some radiation, 
especially for the cases where the kinks are observed at rest, but still kinks exist as localized 
and well defined objects. 



V. CONCLUSIONS 



In conclusion, we have analysed the dynamics of SG kinks in the presence of rapidly 
varying periodic perturbations of different physical nature. We have proposed a rigorous 
analytical approach to derive the "averaged" equation for the slowly varying field component, 
and we have demonstrated that in the main order of the asymptotic procedure the effective 
equation is a renormalized SG equation in the case of the direct driving force or rotating (and 
phase-locked to the external ac driving force) background, and it is a double SG equation for 
the parametric driving force. However, the method itself does allow to find in a rigorous way 
the effective equation for the slowly varying field component in any order of the asymptotic 
expansion in the parameter uo~ l , uj being the frequency of the rapidly varying perturbations 
which has been assumed to be large. 
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Our main purpose was to show which kinds of qualitatively new physical effects may be 
expected in dealing with the renormalized nonlinear dynamics instead of unrenormalized one. 
In particular, we have predicted that the parametric driving force may support oscillations 
of the kink's shape (absent in the SG model) which may be viewed as creation of a shape 
mode of the 2tt— kink characterized by the internal frequency Q so i. For the problem of the 
kink propagation on rotating and oscillating background, we have shown that a periodic ac 
force produces a drift in the kink motion, which may be understood as an effect described 
by an effective dc force to the kink motion in the fremework of an "averaged" nonlinear 
dynamics. 

One of the main conclusions of our analysis and numerical simulations, i.e., that the "av- 
eraged" nonlinear dynamics is drastically modified by rapidly varying (direct or parametric) 
driving force but still may be effectively described by renormalized nonlinear equations, is 
rather general and applicable to many other nonlinear models supporting various kinds of 
solitons. 
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FIGURES 



FIG. 1. The steady-state profile of the 4> (a) and cj) x (b) fields as a function of space in the case 
of a parametric driving force. The parameters are 7 = 0.2, L = 10, uj = 100, and e = 200. 

FIG. 2. The frequency of the internal oscillations of the 2-7r-kink, £l so i, as a function of the 
amplitude of the ac parametric force e for two values of the external force frequency, to = 50 
(squares and crosses) and to = 100 (diamonds and plusses). The results presented by diamonds 
and squares are obtained using the effective double SG equation (^8|) whereas the plusses and 
crosses are the result of direct integration of the parametrically driven SG equation (p3|). Note 
that the agreement between the parametrically driven model (|23| ) and the effective ("averaged") 
model ( |28| ) is better for smaller A = e/uj, when corrections to Eq. (f2g|) from the higher-order terms 
are negligible. 

FIG. 3. The normalized IV curves (for the first Shapiro step) characterizing the SG dynamics 
with the parameters 7 = 0.1, L = 16, e = 12.5, and O = 2.5; and with periodic boundary conditions. 
Shown are the zero-field steps at n = 0, 1, 2, 3, 4, and 5, n being the number of kinks in the system. 
Note that the IV curves clearly cross the zero current axis and the steps are slightly asymmetric 
around the voltage < (ftt >= 2.5; the latter effect is caused by the background oscillations and 
relatively small velocity. 

FIG. 4. The steady-state profile of the 4> (a) and cj> x (b) fields as a function of space over 10 
periods of the external driving force. Parameters are L = 16, 7 = 0.1, e = 12.5, uj = 2.5, and 
/ = 0.25. The value of / is selected at the center of the first Shapiro step and, therefore, the kink 
does not move as follows from the theory because the effective force acting on a kink, / — ^ykco, is 
zero. 

FIG. 5. The same as in Fig. 4 but at / = 0.5. 

FIG. 6. The same as in Fig. 4 but at / = 0. The kink is moving due to the uncompensated 
contribution of dissipative losses —'ykuj. 

FIG. 7. The same as in Fig. 4 but for the case of two kinks introduced by a change of the 
boundary conditions. 

FIG. 8. The same as in Fig. 5 but for the case of two kinks. 
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